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Experiments on tunneling into fractional quantum Hall droplets systematically found tunneling 
exponents smaller than those predicted by the ordinary chiral Luttinger liquid theory. In this note, 
by considering the effects of a smooth boundary, we propose a modification of the theory that 
predicts a reduced exponent, in qualitative agreement with experiments. Such modification consists 
on a self-interaction cubic term, containing higher derivatives of the chiral boson field, that provides 
a non-universal one loop correction to the propagator. 

PACS numbers: 



The study of tunnelling of electrons into both incom- 
pressible and compressible quantum Hall states has been 
the subject of intense research from both the theoretical 
and experimental sides (see fl'] and references therein) . It 
is found that the tunnelling conductivity is non-Ohmic, 
/ oc , with a being a function of the filling fraction v. 
This behavior can be understood within the chiral Lut- 
tinger liquid description of the edge physics advanced by 
Wen in ]^ . Although tunnelling experiments in fractional 
quantum Hall effect (FQHE) systems have shown certain 
degree of agreement with the theoretical predictions ob- 
tained from the chiral Luttinger liquid picture there 
remain certain discrepancies which have been addressed 
by different authors [H, H, 0, 0, H, S Hi ■ In particular, a 
reduction of the order of 10% from the theoretical pre- 
diction for the tunnelling exponent has been observed ex- 
perimentally, as well as the absence of the theoretically 
predicted plateaux structure of this exponent [1, 0, ^. 
These issues have been the subject of intense debate 1| . 

Here we present an alternative derivation of the bound- 
ary effective action by including higher order terms in the 
low energy expansion which are induced by the smooth- 
ness of the boundary. The effects of the edge shape and 
the confining potential on the tunnelling exponents have 
been recently addressed using numerical exact diagonal- 
ization where it was found that tunnelling exponents 
could deviate from the chiral Luttinger liquid universal 
predictions. Our results, obtained analytically starting 
from a microscopic description, are in complete agree- 
ment with these findings. 

More precisely, we start from the microscopic descrip- 
tion of the electron fluid in a strong magnetic field fol- 
lowing and in order to include the effects of a smooth 
boundary, we consider the space dependence of the elec- 
tron density over a boundary region of certain width W . 
This leads generically to a boundary effective theory of 



N coupled chiral bosons, with TV = W/a, a being the 
mean interparticle distance. 

In the simplest case of an almost sharp edge, the effec- 
tive action we obtain is that of a self-interacting chiral 
boson. Using this improved effective action we compute 
the tunnelling exponent which, to one loop order, re- 
ceives a small non-universal negative correction which 
depends on the electron density and the short distance 
cutoff (the " size" of the particles) . By suitably choosing 
the UV cutoff, essentially as the mean distance between 
particles, we show that our results are in agreement with 
the experimentally observed values. Since our theory de- 
scribes gapped bulk excitations, the present approach is 
applicable within the plateaux regions of the quantum 
Hall effect. 

We first derive the boundary effective theory following 
essentially the same steps as in Q but with a smooth 
boundary, the presence of which is shown to be crucial 
e. g. in the computation of tunnelling exponents. Let us 
briefly review the main steps in the derivation: 

The action of a two dimensional system of charged par- 
ticles (charge — e) moving in a constant magnetic field 
Aa — —^e-abx'' and an electric potential V{x°') is 



S = 



p 



dt(^^xf+'^e,,x;x'^-eVix;)^ (1) 



where Xp denote the particle positions, a — 1,2 being the 
space directions and p being a discrete particle index. 

The electric potential will be chosen to be approxi- 
mately constant inside some bounded region M , vanish- 
ing smoothly outside the boundary dM through a region 
£ of width W. This region, in which the electric field 
Ea — daV is non- vanishing, defines the edge of the sam- 
ple. The density po will be chosen, consistently with the 
assumed form of the sample, as approximately constant 
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inside A4 and decreasing within the edge £ to be zero 
outside. 

To see the effect of V{x"') at the edge, we approximate 
it in the neighborhood of each particle by a hnear func- 
tion = —EaX"^, assuming that the length in which 
changes in the electric field are measurable is much larger 
than the cyclotron radius, i.e. r « E/\VE\. Then we 
have an effectively constant electric field around the par- 
ticle and we see that it can be eliminated through the 
boost i° = + where Va — eabEb/B. Then the 
classical motion is given by circular orbits with cyclotron 
frequency eB/m and radii related to the initial velocities 
by r = \xin\Tn / eB , the centers of which are moving with 
velocity —v°' transverse to the electric field, i.e. parallel 
to the edge. 

As a conclusion, for a strong magnetic field the action 
can be written as 

s = Y.^ fdt^abi;4 (2) 
p 

where the new variables Xp are given by the above ex- 
pressions in terms of the position of each particle and 
the electric field in its neighborhood. Note that the ve- 
locity Vp is parallel to dAi, being zero for the particles 
outside £. 

If the mean distance between the particles a = \/Pq 
is small compared to all the external scales involved into 
the problem {i.e. the lengths at which variations on the 
electric field or on the density are measurable), we can 
take the continuum limit. Then, we can replace the vari- 
ables Xp by fields t) where the continuous label 
is fixed by the condition that x°'{y^,t = 0) = y"". Then 
the sum on p can be replaced by an integral 

/p R 
dtd^ypo — eabdti^'i'' (3) 

In the continuum theory, the aforementioned condition 
can be characterized by E/dE >> a and po/dpo >> a, 
where d represents the derivative in some spatial direc- 
tion. Since these derivatives can differ significatively 
when taken in different directions, so does the accuracy 
of the limit. 

To be concrete: at the edge, if the curvature of the 
boundary is small, then we have E/d\\E » a (and the 
same for po), and the limit works in the longitudinal di- 
rection. On the other hand, for a sharp enough edge 
we may have E/d±E ~ a, and we most somehow pay 
attention to the discrete structure on the transverse di- 
rection. We will accomplish this bellow by re-discretizing 
the transverse integration. 

Here it is important to stress that y is the initial po- 
sition of the fluid element in its particular "rest frame" 
in which the electric potential term vanishes. Note that 



the operation of taking the continuum limit introduces 
the function po, which was part of the initial conditions, 
into the action of the system in the form of a measure in 
space. 

The resulting action is invariant under the group of 
measure preserving diffeomorphisms in the y^ plane, i.e. 
changes of coordinates such that y'^ — yi + eijdj^/ po. The 
corresponding conserved generator is 

Q='-^ Jd^y ea6e'^a,i''a,£'^H (4) 

since it is conserved for any function H, its integrand is 
conserved, i.e. 

eate''d,x''d,S:'' = Sabe'^ Oji'd^x^l^^ = 2 (5) 

this relation is a constraint in the values of the fields 
that has to be solved to obtain the independent dynam- 
ical degrees of freedom of the system. 

The next step is to solve the constraint ^ and rewrite 
the action in terms of the dynamical degrees of free- 
dom. Here and in what follows, we will use as parame- 
ters the filling fraction v — pa/eB and the effective non- 
commutativity parameter 6 — l/(27rpo) [9|. 

If we parameterize the field with the help of a new 
field (j> according to 

i'^ ^y'^ + e e'^^dbct) + — e°^e'^ d,(t)djdb(t) + 0{e^) (6) 

then ([5]) is solved up to order 0'^ in the interior of Ai. 
We see in ([5]) that (j) is defined up to the addition of 
an arbitrary function of time — > 4> + g{t)^ and that its 
derivatives account for the deviation of the particles from 
the classical center of their motion. 

At the edge, we take (jB]) as a parametrization of 
which is a solution of the constraint up to subleading 
contributions containing derivatives of 9. 

Replacing this solution in the action ([3]) we obtain 

^ " 8^ ldtd^y''''^''{9ji^9A<l>) + 

+ ^ 9j(e^=4at0a,;09c0)) + O{0^) (7) 

We reintroduce the effect of the electric field by boost- 
ing back y"^ = y"^ — v^t, with the velocity defined above. 
We end up with the Lagrangian description of a fluid, 
where each fluid element is labelled with its initial po- 
sition and it is followed along its motion through the 
plane. Here the velocity is a function of j/', almost con- 
stant within £ and vanishing outside the sample. 

In the special case in which ly, 9 are constants the action 
is a total derivative. The usual procedure is to assume 
this and then impose a boundary to the region in which 
the fluid moves. We choose a different way here, allow- 
ing the parameters to change in space and introducing 
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through them the information about the geometry, with- 
out adding any boundary. A further partial integration 
gives 

S = J dtd^ye'^(d,{iy-^)^d,{dt-v^da)^ + 

Here we see that, wherever the sample is homogeneous, 
the derivatives in front of each term vanish, leaving us 
without any dynamics. On the other hand, the dynam- 
ical degrees of freedom localize at places where there is 
a change in the properties of the material (noncommu- 
tativity and/or filling fraction). This is precisely what 
happens at the edge of the sample, and it is at the core 
of our argument 

To stress this point, let us suppose that the parame- 
ters change as a step function, then their derivatives will 
provide delta functions leading us to a one dimensional 
boundary theory. But we see that our treatment is more 
general in the sense that it includes the possibility of 
smooth edges. 

The action ([8]) is written completely in terms of the 
dynamical degrees of freedom </>. Note that this is not 
a boundary theory because it is defined in the full two 
dimensional space. Nevertheless, the degrees of freedom 
are bounded to the region £ in which there is a change 
in the parameters v and 9 {i.e. in the density), and 
propagate chirally along their level surfaces. 

To proceed further, we define a strip in this setup by 
the condition that the parameters are independent of one 
of the cartesian coordinates = {x,y), say x. Then the 
action becomes 

S = j dtdxdy(dy{iy~'^)(t)dxidt - vdx)<p + 

+ \ dyiOiy-') dy^ {dt - vd,){d,^f) (9) 

To recover the discrete structure in the direction y 
transverse to the edge, we proceed as follows: 

• Cut the space in the y direction in slices at j/„ of 
width Ay ^ a. 

• Replace the integral by a sum over the values of 
the integrand evaluated at each slice, and define 

4>n{x) = (j){x,yn). 

• Replace all the y derivatives by its finite dif- 
ference approximation dy4>(x,yn) — {(pix^yn) ~ 
(/)(x,y„_i))/Ay. 

With this method we obtain an effective theory for the 



fields (j)n{x) living on each slice n, with action 

S = dtdx^^ (^KnCyr{(l)n) + 

+ <(£s(</)„)-f /:i(0„,</.„_i))^ (10) 

where we have defined the free, self interaction and in- 
teraction Lagrangians as 

^J=-i(t>n) = 4'ndx{dt-Vndx)4>n 
^■S{4>n) = 4>n{dt - Vndx){dx(l)nf' 
Cx{(f>n,(f>„-l) = -(l>7i~l{dt - Vndx)idx(l)n)^ (H) 

and the constants k„ and are given by 



1 



1 



1 



2 Ay 



Vn-l 



(12) 



It is important to note that the field 0„ will enter into the 
action only when the constants k„, k'^ are non- vanishing, 
i.e. if there is a change in the properties of the sample 
between the slices n and n + 1 . 

The gauge invariance is now 0„ -f- g{t) adding 

the same g{t) to all the (/)'s. 

Note that when the sample has a sharp edge, i.e. if 
the density is changing in a region whose width is smaller 
that the slicing length W < Ay ~ a, then the whole pro- 
cedure is not applicable and the boundary theory corre- 
sponds to the usual chiral boson theory y. 

Let us suppose that the edge is wider than the slicing 
length, so that we have to keep a finite number of terms 
of the sum in PTI| . n — 1, - ■ ■ ,N. 

We have then a single non-dynamical field -0 = 0_i in 
the interaction term Cx{4>o, 4>-i), all other fields appear- 
ing in the £x(0n,0ri-i) terms being dynamical. The in- 
tegration of this multiplier field will enforce a constraint 
on the field 4>q to which it is coupled. When solved, 
(j)Q = f{x + VQt), and replaced in the action, it implies 
that the interaction term with 0i takes the form 

- ^At'i J dtdxfix + vot)idt - vidx){dx<j}if (13) 

which is different from zero provided that the two ve- 
locities are not equal. Since we have assumed that the 
electric field in the edge zone is approximately constant 
in regions of size a ^ Ay, hence vi « V2 and we can 
discard this term within the present approximation. 

We have then obtained an action very similar to the 
original one but without any Lagrange multiplier 



If ^ 

S = J dtdx^^ (^KnCjr 



+ {Snl - l)/:i(</)„,</'„-l)) (14) 
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where the lower bound of the sum has changed and the 
coefficient — 1) ensures that there is no interaction 
term for the field (j>i . 

Let us consider the case of an almost sharp bound- 
ary, in which the density can be well approximated by 
choosing = 1. In that case the action is 



1 

'8^ 



2 I dtdx [K(/)dx{dt — vdx)4> + 



+ K'(j){dt " vdx){dx(, 



(15) 



where we have dropped the subindexes since we have a 
single field. This is a chiral boson theory with a cubic 
higher derivative self interaction, the latter arising essen- 
tially from the smoothness of the edge. 

To see the effect of this last term in the physical prop- 
erties, we compute the propagator since it is directly re- 
lated to the tunnelling exponent [3]. 

In Fourier space we write for the Feynmann propagator 



G{p, LOp 



27r 



(16) 



^ Kp {iOp — vp) + ie 
while the vertex proportional to 

— k6{p + q + r){ujp — vp) qr (17) 
Then, the one loop corrected propagator reads 



G'ip,u;p) = 



27r 



Kp {Up - vp) - SG{p, UJp) 



(18) 



We consider those Feynmann diagrams whose contribu- 
tion to the corrected propagator is 



SG{p,Up) = ttA^ — (wp - vp)p 



(19) 



where A is a momentum cutoff. Replaced in (|16p this 
gives 



G{p,ujp) 



1 



{k + 5k)p {iUp — vp) 



(20) 



where the non-universal correction to the level is given 
by 



Sk = -TT^A^ 



(21) 



A non vanishing correction to the velocity arises from 
the remaining diagrams but it produces no physical con- 
sequences. 

Since the tunnelling exponent is directly related to the 
level by a w k, the appearance of a correction for the 
latter implies, to first order, a correction to the tunnelling 
exponent. This is our main result. 



To estimate the magnitude of this correction, we need 
to relate the momentum cutoff A to the minimal space 
distance measurable a, which naturally leads to the iden- 
tification A = n/a. Using we then get for this choice 
of the cutoff 



Sk = — - 



(4aVo) 



TT 
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= « -0.196 



(22) 



where we have further identified po — Xjc?. The pre- 
dicted dependence on the density of this non-universal 
correction could in principle be tested experimentally. 

Putting all together this leads to a linear dependance 
of the tunnelling exponent a as a function of \jv which 
is very close to the experimental fit presented in 0] (see 
for example Fig. 3 of this reference). It can be observed 
that the departure of our result from the experimental 
data becomes greater for increasing Ijv^ which could be 
attributed to the need to work with iV > 1 in ((M]) . This 
could in turn be related to the need to consider a wider 
boundary region W > a. The analysis of the higher iV 
case will be presented separately [lO| . 

To summarize, we have constructed a boundary ef- 
fective field theory for a fractional quantum Hall effect 
droplet by including the granularity of the fluid (the in- 
terparticle distance a) and the effects of a smooth bound- 
ary, starting from the microscopic model. By further as- 
suming a space dependent density to define the boundary 
region, we have computed the correction to one loop or- 
der of the tunnelling exponent which is good agreement 
with experimental results for the case of an almost sharp 
boundary. 

Wider boundary regions can be easily treated within 
the present approach, by simply taking higher values of 
N in pl|. 

Other approaches to include a boundary in a Chern- 
Simons description of a quantum Hall droplet have been 
proposed [Uj. In this context, it would be interesting to 



study the connection between the two approaches. 

It should be mentioned that the experimental data in 
Q was reanalyzed in Q in view of the results of [1] , and 
certain degree of agreement with the plateau structure 
was obtained. Here we propose an alternative descrip- 
tion which provides a reasonable agreement with the raw 
experimental data as presented in It would be in- 
teresting to study the higher order theory that we have 
constructed in the present Letter along the lines of the 
approach presented in ref. Q . This could provide a closer 
agreement between theoretical and experimental results. 

It would be also interesting to analyze the conse- 
quences of the higher order corrections induced by the 
smoothness of the boundary that we have obtained here 
in the transition between plateaux states. 

The approach proposed in this note has been formal- 
ized carefully starting from a microscopic description and 
developped further in |12l] . 
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This procedure may have some connection with the so- 
called boundary reconstruction effect, which has been 
studied from both the experimental and analytical point 
of views in 0] . 



